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A B-Spline  Approach  to  Hermite  Subdivision 


U.  Schwanecke  and  B.  Jiittler 


Abstract.  We  present  a new  approach  to  Hermite  subdivision  schemes. 
It  is  based  on  the  observation  that  a sequence  of  second  order  Hermite 
data  define  a unique  interpolating  cubic  C1  spline.  The  B-Spline  form  of 
this  interpolating  spline  leads  to  a stationary  binary  subdivision  scheme 
with  4 different  subdivision  rules  for  the  control  points.  We  construct  a 
generalized  4-point  scheme  which  leads  to  a new  family  of  C 2 Hermite 
subdivision  schemes. 


§1.  Introduction 

Starting  from  an  initial  sequence  {/ifAiezz  of  second  order  Hermite  elements 
(i.e.  vectors  containing  function  values  and  associated  first  derivatives),  a Her- 
mite subdivision  scheme  (cf.  [4, 5,6, 7])  of  order  two  recursively  generates  finer 
sequences  {holies  of  Hermite  elements  associated  with  the  dyadic  points 
pA  = i2-k}i€w  The  refinement  is  based  on  two  rules, 
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where  the  matrix  masks  = {AqA  . . . , A$},  = {J9qA  . . . , B$}  of 

the  scheme  consist  of  real  2x2  matrices  A^,B^  depending  on  the  subdivi- 
sion level  k.  Merrien  [7]  considered  Hermite-type  2-point-schemes  (i.e.  with 
m = 1),  generating  C1  functions.  By  introducing  an  auxiliary  point  subdivi- 
sion scheme,  Dyn  and  Levin  [4,5]  analyzed  stationary  Hermite-interpolatory 
subdivision  schemes  of  arbitrary  order.  Using  this  approach,  Kuijt  [6]  con- 
structed several  C'2  Hermite  interpolatory  subdivision  schemes  of  order  two. 
Kuijt  derived  the  refinement  rules  by  considering  the  polynomials  interpolat- 
ing neighboring  Hermite  elements,  and  sampling  Hermite  data  from  them. 

By  considering  the  interpolating  splines  associated  with  the  Hermite  ele- 
ments, this  paper  introduces  a new  approach  to  Hermite  subdivision.  We  ana- 
lyze the  smoothness  of  the  limit  function,  and  present  a family  of  C2  Hermite 
subdivision  schemes  generalizing  the  4-point  scheme  [3].  This  spline-based 
approach  can  be  generalized  to  Hermite  elements  of  arbitrary  order. 
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§2.  Spline  Subdivision  Schemes 

At  each  subdivision  level  k,  the  Hermite  data  define  a unique  in- 

terpolating cubic  C 1 spline,  having  the  B-spline  representation 

X(fc)(t)  = 5>!fc)ArM(i)  with  knots  t[k),t\k),t (2) 

* 2x  2x 

The  control  points  E 1R  are  associated  with  the  Greville-  abscissas  (see 


e-S-  [8])  tv  = *i 
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sequence.  Control  points  and  Hermite  elements  are  related  by  the  transfer 
mations  (p^),P^+i)T  = H^h\k)  and  h[k)  = (HW)-1(p{£\p%)+1)T , where 
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Clearly,  the  spline  function  can  be  represented  with  respect  to  the  refined 
knot  vector  T^+1b  Knot  insertion  leads  to  the  following  4 refinement  rules 
for  the  B-Spline  control  points: 
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Pii  4^2  i ' 4^2i 
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The  affine  combinations  (4)  describe  a 4-rule  stationary  binary  subdivi- 
sion scheme  for  the  (nonuniformly  parametrized)  B-spline  control  points.  This 
scheme  generalizes  the  splitting  step  of  a binary  uniform  subdivision  scheme. 
The  sequence  of  control  polygons  converges  to  the  C1  limit  function  W°b 
Generalizing  (4)  leads  to  the  notion  of  a spline  subdivision  scheme: 

Definition  1.  A spline  subdivision  scheme  5(a°, a1,  a2, a3)  with  the  coeffi- 
cient masks  ah  = (ap, . . . , o£m+i)>  generating  a sequence  of  cubic  C1  spline 
functions  X^k\t),  is  given  by  the  four  subdivision  rules 

2m+l 

P(4kitk  = E «>2i+i  , h = 0, 1,2, 3 , fc  = 0,1,2, ....  (5) 

3=0 

With  the  help  of  the  transformations  (3),  the  matrix  masks  of  Hermite  sub- 
division schemes  (1)  can  be  transformed  into  the  coefficient  masks  of  spline 
subdivision  scheme  (5),  thus  motivating  the  following  definition. 

Definition  2.  A Hermite  scheme  is  said  to  be  stationary  if  the  matrices  Aj  := 
H(k+V  A(jk)  (ffWf1,  Bj  :=  H(-k+l'>  B{k)  (tfW)-1,  are  constant  for  all  k ( j = 
0, . . . , m).  The  coefficients  ah  of  the  associated  spline  subdivision  scheme  are 
obtained  from 


Consequently,  every  stationary  Hermite  subdivision  scheme  S( is 
equivalent  to  a spline  subdivision  scheme  5(a°,a1,a2,a3).  Note  that  a spline 
subdivision  scheme  can  also  be  seen  as  a special  matrix  subdivision  scheme 
(see  [1])  acting  on  vectors  of  2 consecutive  control  points. 
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§3.  Convergence  Analysis 

In  the  sequel  we  generalize  the  approach  introduced  by  Dyn,  Gregory  and 
Levin  [2]  to  the  spline  subdivision  case.  Consider  a spline  subdivision  scheme 
5(a°,  a1,  a2,  a3)  on  the  finite  domain  [0,n]  € IR.  The  scheme  is  well  defined 
on  this  domain,  for  all  k > 0,  if  the  associated  Hermite  elements  at  stage  k 
are  defined  on  the  set  {t,^1  \i  € Z^},  where 


Zk  — { 0, 1,  - - . , 2kn  + ni}, 


( 2m  - 1 if  Am  0 
| 2m  - 2 if  Am  — 0. 


The  spline  function  (2)  has  the  knots 


T^={t{klik\t{0k). 


,(fc)  ,(fc)  , 

L1  > Ll  ) • * * > 


(*) 
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and  the  control  points  (p|fc))j= 0 2(2*n+n,)+i- 

Consider  an  interval  ljk^  = at  the  fc-th  subdivision  step.  The 

control  points  which  govern  the  future  behavior  of  the  process  in  this  interval 
are  gathered  in  the  vector  pifc  = (p^\  . . . ,p$+ni+1)+1)T. 

The  control  point  vectors  P2i,fc+i , P2i+i,jfc+i  at  the  subdivision  level  k+ 1, 
associated  with  the  subintervals  and  \ are  obtained  from  p,^  by 

two  linear  transformations, 


P2«,fc+i  = Go  p i,k  and  P2*+i,fc+i  = G i p itk  (6) 

with  Go  = 2)  and  Gi  = G(j"^_2),  where  G(V";V)  is  the  matrix 

comprised  of  the  elements  of  G,  at  rows  ii  <■  ■ -<ip  and  columns  j\  <■  ■ ■<  jp. 
These  linear  transformations  are  expressed  as  submatrices  of  the  M x M 
generator  matrix  G,  where  M = 2(ni  + 3).  If  Am  ^ 0,  then  M = 4 (m  + 1), 


and  we  get  the  generator  matrix 
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Otherwise,  if  Am  = 0,  M — 4 m + 2 and  the  generator  matrix  is  as  above  but 
with  the  last  two  rows  and  columns  deleted. 


3.1.  Continuity 

The  following  necessary  condition  is  analogous  to  [2,  Prop.  2.3].  Alternatively 
it  can  be  formulated  using  the  eigenstructure  of  the  masks  of  the  associated 
matrix  subdivision  scheme,  cf.  [1], 
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Proposition  3 (Affine  invariance).  A necessary  condition  for  the  uniform 
convergence  of  the  spline  subdivision  process  to  a continuous  nonzero  limit 
function  on  [0,n],  for  arbitrary  initial  data,  is  that  ak  = 1 for  all 

h = 0,1, 2, 3. 

In  order  to  analyze  the  convergence  to  a continuous  limit  function,  we  exam- 
ine the  difference  scheme  A5(a°,a1,a2,a3)  generating  the  differences  A^  = 
— p[k\  If  the  necessary  condition  of  Proposition  3 is  satisfied,  then  this 
process  can  again  be  described  with  the  help  of  another  generator  matrix 

(8) 

which  is  obtained  using  the  upper  triangular  matrices  Em  — (— Si,j  + 5j+i,j) 
and  {Em)*1  = (—  Shlio1  f>i+h,j),  cf.  [2,  Prop.  3.2].  The  M— 3 differences 
governing  the  future  behavior  of  the  process  in  the  interval  l\k^  are  again 
collected  in  a vector  A,^  = [A^, . . . , A2*)+ni+1j]T.  The  analogues  of  the 
transformation  (6)  are 


^2i,fe+i  = C0Aitk  and  A2i+i,fc+i  = CiA*,*,  (9) 

where  Co  = and  Ci  = ^ii.-M-l)'  Note  that  the  row  and  column 

ranges  of  the  sub-matrices  Co,  Ci  are  different  from  those  in  [2],  as  we  analyze 
a difference  process  with  4 (rather  than  2)  rules.  We  get  (cf.  [2,  Theorem  3.1]) 

Theorem  4.  Let  the  spline  subdivision  process  satisfy  the  necessary  conver- 
gence condition  of  Proposition  3.  Then  the  following  are  equivalent: 

( i ) The  spline  subdivision  process  S'(a0,a1,a2,a3)  converges  uniformly  to  a 
continuous  limit  function  on  [0,  n]  for  arbitrary  initial  data. 

(ii)  The  difference  process  AS(a°,  a1 , a2,  a3)  is  contracting,  i.e.  it  converges 
uniformly  to  zero  on  [0,  n]  for  arbitrary  initial  data. 

(Hi)  There  exists  an  integer  L > 0 and  a real  number  0 < a < 1 such  that 
IlCjj  ■■■CiL  ||oo  < a,  V ij  E {0,1}  and  j = 

In  the  sequel  we  have  to  analyze  other  point  processes  with  four  different 
refinement  rules.  The  continuity  of  the  limit  function  can  then  be  analyzed  in 
an  analogous  way,  where  the  generator  matrix  is  obtained  as  in  (8). 

3.2.  Derivative  process 

In  order  to  investigate  the  differentiability  of  the  limit  function  /,  we  analyze 
the  first  derivative  of  the  cubic  C1  splines  (2).  Clearly,  we  obtain  a sequence  of 
quadratic  C°  splines  with  knots  T;-k\  see  Figure  1.  If  the  necessary  condition 
of  Prop.  3 is  satisfied,  then  the  quadratic  splines  are  generated  by  another 
spline  subdivision  scheme,  again  with  four  different  rules  for  the  control  points. 
This  scheme  will  be  called  the  derivative  scheme  95'(a0,a1,a2,a3). 
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Fig.  1.  Derivative  scheme  and  inscribed  polygon  process. 

Proposition  5.  If  the  derivative  process  dS(a°.  a1,  a2,  a3)  converges  uni- 
formly to  d G C[0,n],  then  the  spline  subdivision  scheme  S( a0,  a1,  a2,  a3) 
converges  uniformly  to  f € C^O,  n],  and  f = d. 

Using  similar  techniques  as  in  Section  3.1,  we  define  control  point  vectors  and 
a generator  matrix  V.  We  omit  the  details,  giving  only  the  main  result: 

Proposition  6.  The  derivative  scheme  dS( a0,  a1,  a2,  a3)  has  the  ( M — 1)  x 
(M  — 1)  generator  matrix 


with  Elf  =diag(l,  2,1,...,  2)EM  and  (Eh)  1 = (EM)  1 diag(l,  §,  1, . . . , |). 

The  continuity  of  the  limit  function  generated  by  the  derivative  scheme 
can  now  easily  be  analyzed  as  in  Section  3.1  by  discussing  the  associated 
difference  scheme  A<9S(a°, a1 , a2,  a3).  This  leads  to  criteria  for  Cl  continuity 
of  Hermite  subdivision  schemes. 

3.3.  Ck  convergence  analysis  via  inscribed  polygons 

In  order  to  examine  higher  order  continuity,  we  inscribe  a polygon  into  the 
quadratic  C°  spline  and  analyze  the  resulting  subdivision  scheme,  called  the 
inscribed  polygon  process  PdS( a0,  a1,  a2,  a3).  More  precisely,  at  the  subdi- 
vision level  k,  we  consider  the  piecewise  linear  function  with  the  vertices 
( t\k+1\ X<'k')(t\k+1'*) ),  see  Figure  1. 

Proposition  7.  The  inscribed  polygon  process  PdS( a0, a1, a2, a3)  of  the 
derivative  scheme  has  the  generator  matrix  P = Lm- i 2?  (1m-i)_1  which  is 
obtained  using  the  ( M — 1)  X ( M — 1)  auxiliary  matrices 
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The  derivative  and  the  inscribed  polygon  processes  are  equivalent: 

Lemma  8.  The  derivative  process  dS( a0,  a1 , a2,  a3)  converges  uniformly  to  a 
continuous  limit  function  f on  [0,  n]  if  and  only  if  the  inscribed  polygon  process 
PdS{ a0, a1, a2,  a3)  converges  uniformly  to  a continuous  function  g £ C[0,n}. 
Moreover,  f = g. 

Proof:  This  equivalence  is  due  to  the  convex  hull  property  of  B-splines,  and 
to  approximation  properties  of  interpolating  quadratic  C°  splines.  □ 

Using  the  inscribed  polygon  process,  we  are  now  able  to  discuss  the  con- 
vergence of  the  spline  subdivision  scheme  to  limit  functions  with  higher  or- 
der differentiability,  by  the  same  analysis  as  in  the  non-Hermite  case  (point 
schemes).  We  simply  have  to  analyze  the  divided  difference  processes 
DvPdS( a0,  a1, a2, a3)  (see  [2,  Theorem  4.2])  of  the  inscribed  polygons,  as 
follows. 

Theorem  9.  If  the  l-th  order  divided  difference  scheme  of  the  inscribed  poly- 
gon process  PdS( a0,  a1,  a2,  a3)  exists  and  converges  uniformly  to  ft  £ (7[0,  n\, 
then  also  the  divided  difference  processes  DvPdS{ a0,  a1,  a2, a3)  exist  and 
converge  uniformly  to  /„  £ Cl~v[0, n]  for  u = 0,1,..., I,  and  = /„. 
Hence,  the  spline  subdivision  scheme  S( a0,  a1,  a2,  a3)  converges  uniformly  to 
g £ C',+1[0,n]  with  = fv. 

For  instance,  in  order  to  prove  that  the  limit  function  generated  by  the  spline 
subdivision  scheme  is  C2,  the  difference  process  ADPdS(a°,  a1 , a2,  a3)  of  the 
divided  difference  scheme  has  to  be  shown  to  be  contractive,  analogously  to 
Section  3.1.  The  divided  difference  scheme  DPdS( a0,  a1,  a2,  a3)  has  the  gen- 
erator matrix  D = 2Em-\  P (Em-\)~1{\ "m-Q-  From  this  matrix  we  get  the 
generator  matrix  C*  = Em-2  D (Em-2)~}  (}  ' M-3)  associated  differ- 

ence scheme  ADPdS(a° , a1,  a2,  a3).  In  order  to  guarantee  a C2  limit  function, 
the  matrix  norms  ||(7t*  • • ■ C*L  ||oo,  V ij  £ {0, 1}  and  j — 1, . . . , L have  to  be 
less  than  1 for  some  L,  where  Cq  = @*  {\  m-s)’  an<^  ^7  = C*(i  M-5)' 

§4.  A Generalized  4-Point  Scheme 

Based  on  a geometric  construction,  Dyn,  Gregory,  and  Levin  [3]  derived  a 
family  of  interpolating  4-point  schemes.  This  family  can  also  be  obtained 
from  an  optimization-based  approach,  as  follows.  Let  the  subdivision  scheme 
generate  a sequence  of  piecewise  linear  functions  Y(k)  with  knots  tt-  and 
control  points  qf\  In  order  to  derive  the  refinement  rules,  we  replace  one 
segment  of  with  two  new  ones  (shown  as  dashed  lines  in  Figure  2,  left), 
subject  to  C°  boundary  conditions.  The  new  vertex  r/^,1 ' is  placed  by  min- 
imizing the  jumps  of  the  first  derivatives  between  new  and  old  polygons.  In 
fact,  minimizing  the  weighted  linear  combination 

= ^[r?+1)(4^i1))-y^+1)(4t+1))]2 
+ 2 w[yifc)(tp))-r|fc+1)(tf))]2  + 2w[yjfc+1)(t^)1)-yifc)(tl(+)1)]2 
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Fig.  2.  Weights  of  objective  functions  in  the  point  (left)  and  spline  (right)  case. 

of  squared  differences  of  derivatives  produces  exactly  the  refinement  rules  of 
the  interpolating  4-point  scheme. 

This  approach  can  be  generalized  to  the  spline  case.  In  order  to  derive  the 
refinement  rules,  we  replace  one  segment  of  X with  two  new  ones  (shown 
with  dashed  control  polygons  in  Figure  2,  right),  subject  to  C 1 boundary  con- 
ditions. The  inner  new  control  points  p^tg^  are  placed  by  minimizing 

the  jumps  of  the  second  derivatives  between  new  and  old  splines.  Minimizing 
the  weighted  linear  combination 

F{pT+\\p %$)  = ^ ix{-+1\4k+1ly)  - #+1)(4m})]2 

+ 8w  [X{_k)(t\K >)— ]2  + [#+1)(4+i)-^+ ' ’(4+1)  l2 

of  the  squared  differences  of  the  second  derivatives  gives  the  refinement  masks 

a0  = (0, 0,  i, 0, 0, 0, 0),  a1  = (0, 0,  ±,  f , 0, 0, 0, 0), 
a2  = (0>  I + — | — 2a>,  | + 3w,  | — 3w,  — \ + 2u>,  | — w,  0), 

a3  = (°>  I - +2w,  | -3w,f  + 3w,-§  - 2w,  § + w,0), 

see  (5).  In  order  to  analyze  C 2 continuity  of  the  limit  function,  we  compute 
the  generator  matrix  of  the  difference  process  ADPdS(a°,  a1,  a2,  a3),  which 
has  to  be  contractive.  Using  the  techniques  of  Section  3.1,  we  estimate  the 
C 2 convergence  range  of  the  parameter  w,  see  Figure  3. 


Fig.  3.  Estimating  the  C 2 convergence  range  of  the  generalized  4 point  scheme. 
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Fig.  4.  Interpolatory  limit  functions  (left)  and  derivatives  (right). 

Two  limit  functions  interpolating  three  given  Hermite  elements  have  been 
drawn  in  Figure  4 (left).  The  functions  have  been  generated  with  the  help 
of  Merrien’s  C 1 scheme  ( a = 0.2,  dashed  curve,  cf.  [4,7]),  and  using  the 
generalized  4-point  scheme  (uj  = 0.015,  solid  curve).  As  can  clearly  be  seen 
from  the  associated  first  derivatives  (right),  the  generalized  4-point  scheme 
produces  a C2  limit  function. 
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